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^ It is sjiown, both anatytically and numerically, that cyclotron resonances can be destroyed in 
dense where (5^is the plas ma frequen cy andd^i* the cyclotron frequency), weakly 

collisional, inhomogeneous plasmas when fl»/^)* 2 rf 1, where v is the collision frequency and r^ 
is the mean Larmor radius. The theor^Ubased upon a model Fokker-Planck equation. It is found 
that the particles make a transition from magnetized to unmagnetized behavior. This is an im¬ 
portant result since it Indicates that the ion- and electron-cyclotron-drift instabilities transform into _ 
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their unmagnetized counterparts, the lower-hybrid-drift instability and the ion acoustic instability, 
respectively. The ion-cyclotron-drift instability (or drift-cyclotron instability) is examined in detail 
and is found to b ecome t he lower-hybrid-drift-instability in the region of maximum growth when 

Q/fij <; Vg/fi ) for T e T t plasmas. The first inequality is required to overcome 

electron viscousaampmgTwhile Ithe secoiul allows the ions to becomfe"unmagnetized.“"‘Applications 
to the equatorial F region of tht ionosphere and the Tandem Mirror Experiment (TMX) are discussed. 
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DESTRUCTION OF CYCLOTRON RESONANCES 
IN WEAKLY COLLISIONAL, INHOMOGENEOUS PLASMAS 

I. INTRODUCTION 

We investigate the effect of like particle collisions (i.e., ion-ion and electron-electron colli¬ 
sion) on short wavelength, electrostatic cyclotron instabilities in weakly collisional, inhomo¬ 
geneous plasmas. That is, for via << 1 and kf rl >> 1 where v is the collision frequency, 
fl is the cyclotron frequency and r L is the mean Larmor radius. We base our analysis on a 
model Fokker-Planck equation originally proposed by Chandresekhar 1 and applied to a homo¬ 
geneous, unmagnetized plasma by Lenard and Bernstein 2 and to a homogeneous, magnetized 
plasma by Dougherty. 3 Dougherty found that cyclotron resonances could be destroyed by colli¬ 
sions if v/(i >1 or (v/a)k? rl> 1. The latter condition is particularly interesting since it 
indicates that cyclotron waves with k} rl »1 may not exist even though via << 1. Physi¬ 
cally this occurs because the particles can diffuse a distance L D — (i//ft) 1/2 r L in one gyro- 
period. If this distance is greater than a wavelength ( i.e., (v/a)kf rl >1) then the cyclotron 
wave cannot maintain its coherence. In fact, it has been shown analytically that the electron- 
Bernstein-mode dispersion equation makes a transition to the ion-acoustic-mode dispersion 
equation ( i.e., the electrons become "unmagnetized" ) when ( vja,)kl rl, > l. 4 We men¬ 
tion that this case has been studied numerically and such a transition was not observed. 3 

We demonstrate in this paper that in weakly collisional, inhomogeneous plasmas, the ion- 
and electron-cyclotron drift instabilities transform into their unmagnetized counterparts, the 
lower-hybrid-drift and the ion acoustic instabilities, respectively. We discuss the ion-cyclotron- 
drift instability (also known as the drift-cyclotron instability) 6 in detail and show it becomes the 
Manuscript Submitted January 31, 1979 
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lower-hybrid-drift instability 7 in the region of maximum growth when (m,/m,) l/2 (w/ft,) > 
vjfl,> m,/m, for T f =sT, plasmas. The first inequality is required to overcome electron viscu- 
ous damping, while the second allows the ions to become "unmagnetized." Applications to the 
equatorial Spread F phenomenon* and the Tandem Mirror Experiment (TMX) 9 are discussed. 

We point out that several Russian authors have also studied the effects of collisions on 
instabilities in inhomogeneous plasmas. Mikhailovskii and Pogutse 10 have derived a general 
dispersion tensor, including electrostatic and electromagnetic perturbations, based on the 
Bhatnagar-Gross-Krook (BGK) collision model." Rukhadze and Silin 12 base their analysis on 
the Landau collision integral and discuss a variety of drift instabilities. However, neither of the 
above studies treat the limit ( v / 0 )k* r£> 1 and hence, do not find a transition to "unmag- 
netized" behavior. 

The structure of the paper is as follows. In the next section we derive a general disper¬ 
sion equation for electrostatic waves in a weakly collisional, inhomogeneous plasma based upon 
a model Fokker-Planck equation. In Section III we apply this theory to the ion-cyclotron-drift 
instability, presenting both analytical and numerical results. Finally, in the last section we dis¬ 
cuss the implications of these results for both space and laboratory plasmas. 

II. THEORY 

A. Physical Configuration and Assumptions 

The physical configuration which we consider is described as follows. The equilibrium 
magnetic field is B - B 0 e x , the density varies only in the x direction and the temperature is 
assumed constant. The density gradient produces a cross-field drift given by l) 0 “ V d i, where 
V d - (v^/20)d In nlbx is the diamagnetic drift velocity. Here, \} h - 2 T/m is the thermal 
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velocity and ft - qBJmc is the cylcotron frequency for a species with temperature T, charge q 
and mass m. We treat only electrostatic oscillations and assume perturbation quantities to vary 
as exp[-/(k - x - a*/)) where k - k ± i y + k u e 2 . We make use of the local approximation which 
requires kL„ » 1 where L„ - (8 In n/Bx)~ l is the scale length of the density inhomogeneity. 
The equilibrium distribution function used in the analysis is 13 


/o - n 0 (x) |^rj ex P 


where we emphasize that this is a self-consistent Vlasov equilibrium for strongly inhomogene¬ 
ous plasmas (i.e., L, > r L where r L is the mean gyroradius) locally, say at x - Xq. 


B. Model Fokker-Planck Equation 

We choose the following model Fokker-Planck equation to describe collisional effects 1 " 3 


¥ + Vj-¥ + 1 (E + i »x B )j-j¥ “ •'-Mfy-t/,)/ + (2) 

8/ ’ bXj m c 1 8v, 8vj I 1 1 m 8vJ 

where n - JVv/, nV, - JVv v,/ and 3nT - JVv m(v-V) 2 / The first term in brackets 
represents friction and is taken to be proportional to the velocity relative to the mean velocity. 
This is a reasonable approximation for thermal particles but is incorrect for "tail* particles. The 
second term represents diffusion and is chosen to be isotropic which requires the Debye length 
to be less than the mean Larmor radius. This model describes small-angle collisions and is par¬ 
ticularly useful in handling like particle collisions. 

Linearizing Eq. (2) we obtain 


-“' a+ £^ 
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«. - $<Pvf\.n 0 Uj - J>v VjA,3n 0 (T\/T 0 ) - (2/v 2 ) J>v(v - U.) 2 /, and vji, - 2 7 ’ 0 /m. 
Here, the subscripts o and 1 denote unperturbed and perturbed quantities, respectively, and 


We can cast Eq. (3) into a form which can be solved using a Green’s function 1 by the fol¬ 
lowing transformation: x—x + V / and v—*v + V where 

v x - . n ~ u 0 \ V y - * U t . 

v 7 + n 2 ' v 2 + n 2 

Making use of this coordinate change we find that Eq. (3) becomes 

i>r +, 0 + |i|i* x "•[- x, '|0- 3x/i -40- <5) 

and h is defined by Eq. (4) with the transformation reflected in /„. We point out that the above 
coordinate transformation is sign dependent (through the diamagnetic drift velocity (/„ ) and 
cannot be made simultaneously for both species. However, this does ot create a serious prob¬ 
lem since we only consider collisional effects on a single species. Moreover, since we are 
interested in the situation where v/ft« 1 and kj r£ » 1 we note that (1) the old reference 
frame is a good approximation to the new one and (2) the final term in Eq. (4) can be 
neglected. 4 Thus, within the context of these assumptions, we consider 

Hr * v 4 + 1 i (r * ■•), - "'10 - Wi - 40 - i 00 <6) 

to describe electrostatic oscillations in a weakly collisional, inhomogeneous plasma where /„ is 
defined by Eq. (1). The important physics contained in Eq. (6) is the modification of the 
unperturbed orbits to include a resistive term and the diffusion of perturbations as they are con¬ 
nected along an unperturbed orbit. 

C. Derivation of the Dispersion Equation 


The dispersion equation for electrostatic perturbations is 


D(«. k) - 1 + Zx. -0 


(7) 
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where x„ is the susceptability of species a is defined by x„<t> i “ - (4ire„M 2 )/,„,. Thus, we 
only need to know the perturbed density in order to determine the linear behavior of electros¬ 
tatic waves. In order to compute the perturbed density we follow the method of Dougherty and 
employ a Green’s function. 3 We only present the results here and refer the interested reader to 
Ref. 3 for the mathematical details. In the absence of collisions, the Green’s function is a 8 
function and perturbations are simply convected along the unperturbed orbit. When collisions 
are included, the perturbations also spread out due to the diffusion as well as being convected 
along the unperturbed orbit (which now includes a resistive term). Thus, the Green’s function 
becomes a Gaussian distribution in phase space. 


Using notation similar to Dougherty, we find 

", - ~T «'* 


HvL- 

/ - drtxpli&r - <Mt) - V(t)), 

4>(t) - y kl rl [cos# + vr - exp( - Pt)cos(t - 0)1 

t *,? r? 

+ ——— 1 — exp( - *>r)] 

2 i 

t y o 

+ ; - ^ - [exp( - i»r)sinT + ?(1 - exp( - Pt)cost) 1, 


——- _ -+ <r K p x exp[ - (i> + ,X)t] 

2/o» Itf 


4R 

P x 1 - exp[ - (?-f /X)rl 
v + iX 

expl - (» + /X)t!}, 


a» -- 2 -^, (Mr* - «f“), 


( 10 ) 

( 11 ) 


(12) 
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9 - 2tan~‘ P.P - v/(l,<L - w/ftandej 1 - (dlnn/dx^),., . Note that we use a polarized coordi¬ 
nate system Or 1 ,* 0 .* -1 ) where x' - (l/yft)(x + iy),x° - z and jr' 1 - (1 /a/D(x - iy). This 
transformation is unitary but not orthogonal so we must distinquish contravariant and covariant 
vectors. 


Making use of Eqs. (8) - (12), we arrive at 

X “ ( 1 + '{ ft U ° jX </Tex Pf , ‘ iT “ < P(’’)l| (13) 

where - 4 mnq l !m. We emphasize that this expression for \ is only valid for i//ft «1 and 
kl rl» 1. Note that in the limit £,„—«> ( i.e., U o —0 ), Eq. (13) reduces to a result 
equivalent to that derived by Allan and Sanderson 4 and in the collisionless limit (P — o) 
becomes 


I m-k ± U, 


which is a standard result. 15 


J* dr exp I itor + y/tj 2 r t 2 (cost — 1) — / -^yp-sin r (14) 


I- j *n rli 2 


III. APPLICATION: ION-CYCLOTRON-DRIFT INSTABILITY 

The ion-cyclotron-drift (ICD) instability (also known as the drift-cyclotron instability) 5 
has received considerable attention in the past IS years. It is thought to play an important role 
as an anomalous transport mechanism in a variety of "collisionless" plasmas (e.g., post¬ 
implosion stage of a theta pinch, 15 solenoidal section of TMX, 17 the polar cusp, 1 * the iono¬ 
sphere, 19 the Earth’s magnetotail 20 ). The instability can be excited when LJr L < (m,/m,) l/2 // 
and is characterized by w, — 1(1, and y — I(l ,(m,/m,) l/4 where <o - w r + iy and / is the ion 
cyclotron harmonic number. Maximum growth occurs for It • B - 0 (i.e., - 0) with 
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kr Le ~ (T t /T ,) u2 in dense plasmas (i.e., u 2 , » ft 2 ). Note that to satisfy the condition 
(v ,,/ft i)k*rfr > 1 at maximum growth, we require 



For a hydrogen plasma, Eq. (15) implies that if (i/,,/ft,) > 5 x 10 -4 , ion-ion collisions will 
destroy the cyclotron resonances and hence, the ICD instability. However, rather than drive 
the plasma to stability we demonstrate, both analytically and numerically, the transformation of 
the ICD instability to the lower-hybrid-drift (LHD) instability in the presence of weak ion-ion 
collisions. 

The dispersion equation under consideration is 

Dim, k) - l + Xi + X'~ 0 (16) 

where 

*' “ j 1 + *1 “ n ^ |X ^expOwr/ft, - 4 >,(t))J, (17) 

2ti> 2 f I ( y II 

Xe - ll —exp(6 f )/„(6 ( ,)ll—— , (18) 

and 

♦,(r) - jfcj 2 rfc Jcoso- + i>, r—exp(—i>, t)cos(t—8) J (19) 

+ i ~~’~ jexp(-P,T>sinr + J, (1 - exp( - P,t)cost)J, 

V, - v„/ft,. b c - j- k : ri„ g» 2 „ - 4 nnellm a ,vl - 2 TJm a , V da - (v a 2 /2ftj(d In «/dx),_, 0 
ft,, - e a B(Jm a c, I, is the modified Bessel function of order /. The simple form chosen for the 
electron susceptibility is valid since ft, < at « ft, and it is assumed that v e J(i t « 1, 
p„/ft, « 1 and v„k J r £ « ot (i.e., the electrons are collisionless). The final condition is 
required to avoid collisional damping due to electron viscosity 21 and will be discussed further in 
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the final section. We emphasize that Eq. (17) is valid only in the regime 
vj(l, « 1 and klrl » 1. 


A. Analytical Results 

i. v(t?rl) « 1 

In the limit PUc 2 r£,) « 1, v « 1 and tfrfc » 1 collisional effects can be ignored and 
we find that 



2 [. 1 1 ^ « - kV d , | 

*■" l 1 "7? 122 

where (VwAr t ,) -1 is the large argument expansion of /_ |y X : r/ J exp | - y * 3 r£ J. The ioi 


cyclotron resonances are apparent in Eq. (22). 


ii. PUV,) > 1 


In the limit vik 1 ^,) > 1, v « 1 and k 2 r£, » 1 we can again approximate the electron 
susceptibility by Eq. (21). However, the ion susceptibility is changed substantially from Eq. 
(20). We consider Eqs. (17) and (19) and point out that for k? r£,» 1 the dominant contri¬ 
bution to the integrand occurs in the vicinity of T=2irm(m-0, 1, 2,..). We let 


2 nm + i and obtain 
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/„ dr expOwr - ®,(t)) = dt exp j/ j- - - - —| t - y * 2 r t 2 / 2 j (23) 

+ Z m<tt exp [' |~ n *-- -| t - y tfrlt 1 J exp[-2rr* 2 r£«?] exp|+2irwi 

where 8o, 8 m « 1. Clearly when H&rfr) > 1, the first term in Eq. (23) is dominant and we 

find (letting 8 q — «>) that 

\2 Wp 

Xi ~ (1 + ( 24 ) 

where {, - (w - kV dl )/k\, and we have made use of the relation 

Z({) - iS/ir exp(-C 2 )(l + er/(/£)l. (25) 

Thus, the dispersion equation becomes 

D(w. *) - 1 + [1 + {,Z(£,)1 + ^ j 1 - J 1 - -y^| exp(-6,)/ 0 (6,)j (26) 

which describes the lower-hybrid-drift instability. 7 Physically, the instability is excited since the 
ions can now move across the field lines and be in resonance with the drift wave. We mention 
that this analysis parallels that of Allan and Sanderson who considered electron Bernstein 
modes in a homogeneous plasma. 22 Finally, note from Eq. (23) that the ICD instability will 
also transform into the LHD instability in the strong drift velocity regime (i.e. y > (1,). 

B. Numerical Analysis 

We now solve Eqs. (15)-(18) numerically. As a relevant example, we consider the equa¬ 
torial F region of the ionosphere where the dominant ionic component is O + so that 
m, -16 m p . Recent experimental observations indicate intense VHF and UHF radar back- 
scatter during equatorial spread F resulting from density irregularities of 1 m and 36 cm, respec¬ 
tively. 19 These irregularities correspond to wavelengths such that k 2 rfr » 1. Since there is 
also evidence of density inhomogeneities with scale lengths L„/r Ll < (m,/m,) U2 , 23 it has been 
suggested that the ICD or LHD instabilities are responsible for the irregularities. Also, because 
v,t/Sl, < 10~ 2 ion-ion collisions can play a significant role and a numerical study is required. 

9 
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In Fig. 1 we show y/ft, vs. kr u and kr L , for typical ionospheric parameters: T t - T„ 
«y/ft, — 10.0, F rf( /v, — 0.037 which corresponds to L„ ~ 13 r L/ (~ 75 m) and 1 ',,/ft, — 0.0, 
10 -6 , 10 -5 , 10 -4 . We comment that wave growth also occurs for kr^ > 1.5 but has not been 
plotted. We first indicate the wave growth spectrum for the LHD instability (—) where we 
have used Eq. (25). The spectrum is broadband with maximum growth at kr u = 1.2 with 
w, ~ 3.0ft, and y ~ 0.08 ft,. Growth occurs for kr Lf < 0.15 but is very weak (i.e., 

y < 10 _3 ft,). For •',,/fl, - 0.0 (-—) we see that the first three cyclotron harmonics of 

the ICD instability are excited with maximum growth occurring at kr u ~ 1.0 with u r = 3.0 ft, 

and y ~ 0.15 ft,. As t^/ft, is increased to 10 -6 (-) we already note a change in the 

wave spectrum. Although the harmonic structure is maintained, the maximum growth rate 
decreases 6% and growth is extended beyond kr u — 1.45 for the 3rd harmonic. For 

*,,/ft, - 10~ 5 (.) a more dramatic change occurs in the wave spectrum as the harmonic 

structure becomes distorted. The wave spectrum collapses to the LHD curve for kr u ~ 1.4 
( kr Lt = 210) which corresponds to (i» // /ft / )fc 2 r^ — 0.5. Moreover the maximum growth rate 
decreases by 40% to y ~ 0.09 ft,. The first two harmonics are still evident although their 
growth rates have decreased somewhat. When 1 ',,/ft, - 1(T 4 (xxxxx) the wave spectrum falls 
on the LHD curve and only the slighest hint of harmonic structure remains at / - 1 and / » 2. 

Figure 2 plots ,'^/ft, vs. VJv, for T t - T., oy/ft, - 10.0 and v,,/ft, - 0. 10" 4 , 10"\ 
10 -4 . Here y M is the growth rate maximized with respect to k. Again we indicate the results 
of the LHD instability (-) as a reference. Setting i',,/ft, - 0 we recover the ICD instability 
which displays substantially stronger growth than the LHD instability for weak drifts (i.e., 
v d , < 0.10 V,). We can distinguish the excitation of the first 8 harmonics as VJv t is increased. 
Also, note the transition to the LHD instability for V dt >0.11 ft, which occurs when 
y — 0.75 ft, as indicated earlier. For 1 ',,/ft, - 10"‘ (-) we note that the growth rates 
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are generally reduced for K* < 0.08S v, although 8 harmonics are still apparent. Increasing 

*„/n, to 10~ s ( . ) causes the wave spectrum to collapse to the LHD spectrum with only the 

first S harmonics visible. Finally, for v„/0, - 10~ 4 the curve lies on the LHD plot. 

Thus, we have shown numerically that the ICD instability makes a transition to the LHD 
instability for > 1. Since growth occurs for k 1 ^, » 1, we find that only a very 

small amount of collisions j-~p > 10~ 5 j is needed to cause this transition in the region of 
maximum growth. For typical ionospheric conditions, 10* cm" J < n < 10 4 cm* 1 , 7; — 0.1 eV 
and B — 0.3 C we find that 10~ 4 < -j-p- < 10"' so that the LHD instability should dominate 

when instability occurs. However, it is found that the condition v„ k 1 rl, < m is not Mtisfled 
when n > 10 4 cm -1 for moderate density gradients > 8.0 r L ,) and it is expected that the 
LHD instability will be stabilized by the viscuous damping of the electrons for this situation. 
Hence, instability is likely to occur infrequently and will be restricted to low density regions of 
the F region, such as the interior of ionspheric bubbles during equatorial Spread F**. We defer 
a more complete discussion of this instability and its application to Spread F to a subsequent 
report. 

IV. DISCUSSION 

We have demonstrated, both analytically and numerically, that cyclotron resonances can 
be destroyed in dense («, > ft), weakly collisional, inhomogeneous plasmas when 
(a/f \)k*rl £ |. The physical mechanism for this phenomenon is that particles can diffuse a 
distance L D — (*/ ll)' n r L in one gyroperiod due to collisions. If this distance exceeds the 
wavelength (i.e., (a/O )k 2 r[ > I) then the wave can no longer maintain its coherence. This 
result is in accordance with the work of Dougherty 1 and Allan and Sanderson 4 who considered 
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only homogeneous plumes. It is an important result since it indicates that the ion and electron 
drift-cyclotron instabilities transform into their unmagnetized counterparts, the lower-hybrid- 
drift (LHD) instability and the ion acoustic (IA) instability, respectively, for weakly collisional 
plumu. We specifically discussed the ICD instability in detail and found the transition to 
occur in the region of maximum growth for > ( m t /m ,). However, we have also 

pointed out that electron viscous damping can stabilize the LHD instability when 
v„ k 2 rl, > w since electron viscosity provides a dissipative mechanism for the wave. Recog¬ 
nizing that »}\ we find that (m t /m,) ul «/0, > v„/n, is required for instabil¬ 

ity to occur in the region of maximum growth. Electron-electron collisions then have the effect 
of placing a threshold condition on the diamagnetic drift velocity necessary to excite the LHD 
instability since «— kV#. Hence, 'collisionless” plumas on cyclotron time scales, such u the F 
region in the ionosphere, may in fact be strongly affected by weak collisions as we have shown. 

In collisionless plasmas both the ion and electron cyclotron drift instabilities require a crit¬ 
ical density gradient scale length to be excited. This occurs because the instabilities are pro¬ 
duced by a coupling of a cyclotron wave «i - ft and a drift wave <u 2 ” kV d , and a critical drift 
velocity is therefore required to satisfy the matching condition <o,= <u 2 . However, when 
(v/H)k 2 r£ > 1 these instabilities transform into the LHD and IA instabilities which are driven 
via wave-reticle resonances and have different turn-on criteria. In the case of the LHD insta¬ 
bility a critical drift velocity is required to overcome electron viscuous damping as discussed 
above. On the other hand, the critical drift velocity for the IA instability is sensitive to T,/T t 
and may or may not be excited depending on the actual plasma conditions. 

The difference in the nonlinear behavior between the ICD and LHD instabilities, and the 
ECD and IA instabilities is probably a more important consequence of the theory since entirely 
different stabilization mechanisms are operable. The ECD instability hu been studied in detail 
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and it has been found that in the collisionless limit, resonance broadening due to turbulence 
effectively "unmagnetizes* the electrons and a transition is made to the IA instability. 26 This is 
the same result we have found for (p w /n,)**/£ > 1. A variety of theories have been pro¬ 
posed for the nonlinear saturation of the IA instability and they will not be discussed here. 

Recent nonlinear theories have been presented for the stabilization of the ICD instability 
and rely upon a nonlinear frequency shift 27 or low-frequency (<u « (l,) density fluctuations. 21 
On the other hand, the suggested saturation mechanisms for the LHD instability are quasilinear 
stabilization, 24 ion trapping 20 and electron resonance broadening. 11 We comment that 
Hasegawa’s theory 21 concerning stabilization via low-frequency density fluctuations is also appli¬ 
cable to the LHD instability. We will not discuss the details of these stabilization mechanisms 
but simply emphasize that different saturation energies are anticipated for the ICD and LHD 
instabilities. Hence, the role of weak ion-ion collisions can substantially alter the nonlinear 
plasma response due to different anomalous transport properties. 

Although we have only discussed an application of this theory to the ionosphere, we wish 
to point out its relevance to the Tandem Mirror Experiment (TMX). 9 . TMX will consist of a 
straight solenoidal plasma which is confined by conventional mirror end cells. The solenoidal 
section will be similar to a theta pinch and the ions are expected to have a roughly Maxwellian 
distribution. Since density gradients are expected with scale lengths such that 
Lf/r Ll < it is anticipated that the ICD instability is likely to occur. Using expected 

operating parameters for the solenoidal section (7J - 80 eV, n = 10 11 cm -1 , B - 500 G), we 
And that — 4 x 10~ 3 . Thus, we conclude that drift waves with kr u > 15 will probably 
be excited by the LHD instability and the nonlinear dynamics in this regime can be studied 
using unmagnetized ions. Of course, finite p and temperature gradient effects are also expected 
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to pl«y a role in TMX 17 - ,2 - ,:1 which are not considered in this paper. However, in principle 
these effects can be included in a straightforward manner and are not expected to alter this con¬ 
clusion. 

Finally, we point out that the Green’s function method used in this work is not the most 
direct technique, although it lends itself to a simple physical interpretation. Recently, Catto has 
computed the perturbed density in a collisional, inhomogeneous plasma using a velocity 
transform method 74 which is a simpler analysis. For a more complex situation (e.q., inclusion 
of temperature gradients, V B and electromagnetic effects ) this procedure would be preferable 
to the one outlined in this work. 
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